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Abstract 

In this work we study the solvabihty of the initial boundary value problems, which model a 
quasi-static nonlinear behavior of ferroelectric materials. Similar to the metal plasticity the energy 
functional of a ferroelectric material can be additively decomposed into reversible and remanent 
parts. The remanent part associated with the remanent state of the material is assumed to be a 
convex non-quadratic function / of internal variables. In this work we introduce the notion of the 
measure- valued solutions for the ferroelectric models and show their existence in the rate-dependent 
case assuming the coercivity of the function /. Regularizing the energy functional by a quadratic 
positive definite term, which can be viewed as hardening, we show the existence of measure-valued 
solutions for the rate- independent and rate-dependent problems avoiding the coercivity assumption 
on /. 

1 Introduction and setting of the problem 

Due to the ability of ferroelectric materials to transform a mechanical action into an electrical impulse 
and vice versa they are being used in a broad range of modern engineering devices as actuators and sen- 
sors. Recent technological developments enabled the reduction of the production costs for ferroelectric 
ceramics and thereby increased the interest to use them in the novel implementations. Demand for the 
reliable mathematical models, which on the one hand are capable to describe a complicated nonlinear 
electromechanical behavior of ferroelectric devices in order to optimize their design and predict failure 
processes and on the other hand are simple enough for numerical implementations, caused a rapid 
progress in this field in the last years. A brief review of recent advances in modeling of ferroelectric 
material behavior can be found in [7J. In the present work we study the solvability of the nonlinear 
initial boundary value problems associated with phenomenological constitutive models of ferroelectrics 
[31 m El [HI [m mi ■ similar to models in the metal plasticity the type of ferroelectric models consid- 
ered here is formulated within a thermodynamic framework by using the standard material relevant 
description method of an energy function and a flow rule. In contrast to the micro-electromechanical 
models, which contain a large number of internal variables standing for the distribution and the volume 
interaction of ferroelectric domains, the main goal of the phenomenological models mentioned above 
is to improve the speed and the robustness of numerical implementations by keeping the number of 
internal variables as small as possible. The models presented in [21 [H [51 [HI \IM [HI ^-nd studied in this 
work use as internal variables only the remanent strain and the remanent polarization. 



Setting of the problem. The model equations are formulated as follows. Let 51 C M'^ be an open 
bounded set with the C^-boundary Sfi and denote the set of symmetric (3 x 3)-matrices. Unknown 
are the displacement field u{t,x) € R^, the Cauchy stress tensor a{t,x) e S^, the remanent strain 
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tensor r{t,x) E S^, the electric potential (p{t,x) E M, the vector of electric displacements D{t,x) E M'^ 
and the vector of remanent polarization P(t, x) E MP in a material point x at time t. The symbols 

e{u{t,x)) = ^{\/,u{t,x) + {^,u{t,x)f ) E 

and 

denote the linearized strain tensor and the electric field vector, respectively (e and E for short). The 
fundamental assumption of the models under consideration is that the strain tensor e and the vector 
of electric displacements D can be additively decomposed into reversible and irreversible parts, i.e. 

e = (e - r) + r, D = {D - P) + P. 

In this case e — r and D ~ P are reversible and r and P are irreversible parts of e and D, respectively. 
For {t,x) E fir '■= (0,T) x 51 the unknown functions satisfy the following system of equations 

— div (7 = 6, (1.1a) 

div D^q, (1.1b) 

a = C{e-r)- e^E, (1.1c) 

D^e{e-r)+e E + P, (l.ld) 

''^^dj(:~i^-c(:)) (Lie) 



Ptj -^yyE-fpj yp 

completed by the initial conditions 

r(0, x) = r"{x), P(0, x) = P°{x), xEft (l.lf) 
and the homogeneous Dirichlet boundary conditions 

u{t,x) = 0, (i){t,x) = 0, {t,x) E [o,r) X dn. (i.ig) 

The equations p.lap and (|l.lb[) are the equilibrium equation and the Gauss equation in a quasi-static 
case, respectively. Here, the function b{t, x) E K."^ denotes a given body force and q{t, a;) G M is a given 
density of free charge carriers. The functions g, f : xR'^ ^ M in (|l.le[) denote constitutive functions, 
the form of which are usually determined by experiments. Based on the thermodynamical considerations 
we give in the next two paragraphs the precise conditions, which g and / should satisfy, and discuss the 
equation (|l.lep . The mapping £ : S*^ x — )• S'"^ x M"^ in equation (jl.lep is linear symmetric and positive 
semi-definite and stands for the hardening effects. This mapping is not contained in the engineering 
models considered here and is introduced because of mathematical reasons, which are discussed in the 
last two paragraphs of the introduction. An overview of the previous results concerning the existence 
theory for the ferroelectric models and the structure of the present work can be found in the last 
paragraph of this section as well. 

Due to the additive splitting of the strain tensor e and the vector of electric displacements D into the 
reversible and irreversible parts, the constitutive relations (jl.lep . (jl.ldp can be equivalently rewritten 
as follows 

a= (C + eTe-ie)(£-r)-eTg-i(^_p) 

E= -e-^e{e-r)+e-^{D- P), ^ ' 

that implies that the reversible parts of e and D satisfy the constitutive equations of linear piezoelec- 
tricity. Here the mappings C : ^ S^, e : R"^ — ?> R"^, e : S''^ — >■ R"^ are material dependent elastic, 
dielectric and piezoelectric tensors, respectively. In the engineering literature [21 IH [SI IHl [121 Ull the en- 
tries of the constitutive tensors C, e and e often depend on the internal variables r and P. For example, 
in [8] the tensor e has the following form 

|P| 1 

^kij = —{essnknirij + e^nUkUij + -ei^iniajk + rijaik)), (1.3) 
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where n — jp^, a-ij — Sij — niUj, Cij are constants and Pg is the remanent polarization saturation 
constant. However, because of the difficulties arising in the mathematical treatment of the problem 
()f .ip . in the present work we suppose that the tensors C, e and e are independent of the internal variables 
r and P. Our approach to the derivation of the existence of the solutions for (jl.ip relies heavily on 
the L^-existence theory for elliptic systems with 2 < p < oo. In order to apply such a theory to our 
purposes we have to require that the entries of the tensors C, e and e are continuous functions of a; € f2. 
But since it is expected that the functions r, P belong only to LP{U) for some 2 < p < oo one can not 
guarantee that the mappings C, e and e possess this regularity. Therefore, we suppose that the tensors 
C, e and e are independent of P and r and continuous functions of a; G ri. Additionally, according to 
the engineering models considered here we assume that the mappings C, e and e are linear and bounded 
and that C and e are symmetric and positive definite uniformly with respect to x Q. 

The method presented in this work can be easily generalized to the case of nonhomogeneous Dirichlet, 
Neumann or mixed boundary conditions. 

Thermo dynamical considerations and choices of the function g. A general form of the energy 
function corresponding to the models considered here can be derived by using the constitutive relations 
p.2p and the Clausius-Duhem inequality. Although the different types of thermodynamic potentials are 
used in the literature (for example, the Helmholtz free energy function in [51 [7] or the enthalpy function 
in [12J) it is typical in modeling of the nonlinear behavior of ferroelectric materials to derive model 
equations by means of the Helmholtz free energy in the form '5 = ^'(e, r, P). The main requirement 
is that the function ^! satisfies the Clausius-Duhem inequality 

< ere + (cr - ^,)e + {E - fz,)!) - - *p-P. (1-4) 

The arguments in the thermodynamics of irreversible processes yield that the equations 

cr = ^,{e,D,r,P) and E = ^D{e,D,r,P) (1.5) 

hold. The Clausius-Duhem inequality can be then reduced to the following inequality 

<-^rr--9pP. (1.6) 

Integrating the relations ()1.5p and using (|1.2p we conclude that the free energy function can be repre- 
sented in the form 

4'(e,D,r,P) = *,ev + /(r,P), (1.7) 

where ^-^v = | {{C + e){s ~ r), e - r) - {e^e-^{D ~ P),e ~ r) + i {e-^D ~ P), D - P) is the 

reversible part of the energy. The function / corresponds to the remanent state of the material under 
consideration and is given by 

/(r,P) = /(r,P) + i|/:i/2(^^p)T|2_ 

The authors of the engineering models [Sj El HI [121 [Ml make different assumptions concerning the form 
of the function /. Their choices are usually based on the experimental results. Several examples of the 
function / are given below. The quadratic term with the linear positive semi-definite operator C is not 
contained in the models considered here. It can be regarded as a hardening term and the reason of its 
introduction is discussed in the next two paragraphs. 

Since the entries of the given tensors C, e and e are assumed to be continuous functions and inde- 
pendent of r and P, using the expression for ^ we rewrite the Clausius-Duhem inequality (|1.6p as 
follows 

< (^ - fr)r + {E~ fp)P. (1.8) 

The second law of thermodynamics (|1.8I) restricts the choice of the function g in the equation (|l.lep . 
The inequality p.Sp holds if 5 is a proper convex function. Additionally, we suppose that the function 
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g is lower semi-continuous. In most models in [3J [HI HI [HI III] the function g is chosen as an indicator 
function of some bounded, closed and convex set K C x M.'^ with E K, namely, 




(1.9) 

This choice of the function g corresponds to a rate-independent process. Rate-dependent effects such 
as time-dependent relaxation of ferroelectric polycrystals have been observed experimentally as well. 
To describe the rate-dependent behavior of a ferroelectric material in [6J the function g is chosen in the 
form of a polynomial. In the rate-dependent case we require that the function g satisfies the following 
two-sided estimate with ci, C3 > and 02,04 > 

ci\v\P - C2 < g{v) < oslvl'' + C4, (1.10) 

which holds for any v & S"^ xR^. The condition p.lOp implies that 

g*{v)>di\v\p' -d2, (1.11) 

for di > 0, £^2 > and any v € x M.^, where g* is the Legendre-Fenchel conjugate of g (see 
Appendix |X] for basics on convex analysis) . Throughout the whole work we assume that the number p 
satisfies 2 < p < 00 with p* such that 1/p + 1/p* = 1. 

Possible choices of the function /. In most models in the engineering literature the remanent 
part of the energy / : S*"^ x M'' — > M is given by a convex function whose domain dom(/) is a convex 
(possibly unbounded) open subset of xM.^. In [T21[I11[S] it is assumed that the function f — f{r, P) 
depends only on P. In particular, in [TT, 11] / has the following form 

;(,p) = ;(p)=|f ((l + ^)Ml + ^) + (l-^)ln(l-^)), ma)\<P., 
l+oo, ma)\>Ps, 

where a S is a given direction with ||a|| — 1 and Pg is a saturation constant, and in ^ the function 
/ is of the form 

/(.P,=/(P,J-^?('°('-W> + W)- (1,13) 

[+<X, \P\>Ps- 

In Theorem l2.5l we suppose that if the function / depends only on P, then it has to satisfy the following 
coercivity condition 

/(r,P) = f{P) > ai|F|2 - a2, ai > 0, 02 > 0. (1.14) 

The coercivity condition (I1.14p is satisfied by the function / given by (|1.13p , but not by the function / 
in (|1.12p . Thus, the result of Theorem 12.51 can not be applied to the function / defined by (|1.12l) . 

In [6] and [3] the function / depends on both internal variables r and P and satisfies the following 
coercivity condition 

f{z)>bi\z\P~b2, &i>0, &2>0, z = (r,P). (1.15) 

The present work is especially focused on the rate-dependent processes with the function g satisfying 
the polynomial growth condition (jl.lOp . Under the condition (|1.10p we prove the existence of the 
measure-valued solutions in the sense of Definition 12.21 for the model (|l.ll) without the regularizing 
term, i.e. with C = (see Theorem 12. 5p . However, in this case we have to assume that the function / 
satisfies one of the coercivity conditions given above, i.e. either (|1.14p if / depends only on P or (|1.15p 
if / depends on both variables r and P. If the linear mapping C in (jl.lep is positive definite, then we 
are able to prove the existence of measure- valued solutions for the problem (jl.ip in the rate-dependent 
case without assuming the coercivity of the function / (see Theorem 12. 6p as well. 
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To prove the existence of measure- valued solutions for the rate-independent case (Theorem 12. 6[) 
the linear symmetric and positive definite mapping £ : 5'' x M^^ — > S''' x M'^ in the equation p.lep 
is introduced. As we mentioned in the previous paragraphs this mapping is not contained in the 
engineering models. The additional quadratic hardening term (£z, z) with z = (p) in the energy 
function 'i'{e, D,r, P) regularizes the model p.ip with g given by (|1.9p and the existence of measure- 
valued solutions in the sense of Definition 12.21 can be obtained. The well-posedness of the problem 
in the rate-independent case without the regularizing term C is an open problem at the moment. 

Previous results and structure of the present work. The first existence result for the nonlinear 
ferroelectric models in the rate-independent case is obtained in [lOj via the energetic approach. However, 
in order to use the compactness argument the authors of [TU] regularize the energy function by the 
additive quadratic gradient term of internal variables (£Vz, Vz) with z = (r, P) and positive definite 
£, which is not present in the energy function (|1.7p . For such a modification of the model the authors 
of [TO] prove the existence of strong solutions. Hereby the tensors C, e,e are allowed to depend on 
the internal variables. For the free energy regularized by the term {Cz,z), which can be regarded 
as the hardening, we mention the following existence results O [9]. In these works the tensors C, e, e 
are independent of the internal variables. The derivations of these results require that the Nemytskii 
operator F : L'^{fl) — > R, generated by the function / : 5'^ x — >• R, is Frechet differentiable in 
L^(ri). The last requirement is satisfied if and only if the function V/ is affine. In [Sj^Qj the existence 
and uniqueness of the strong solution is shown in the rate-independent case under the assumption 
F e C2'"P(L2(f7,S'3 x R3)) and F e C^{L'^{n,S^ x R^)), respectively. In the rate-dependent case it 
is believed that there are no mathematical results concerning the existence of solutions. In the present 
work we show the existence of measure- valued solutions of the rate-dependent problem when 
C = 0, i.e. the energy function in (jl.7p does not contain regularizing terms. In Section [5] we introduce 
and motivate the notion of measure-valued solutions for the ferroelectric model formed by equations 
(|l.ip as well as formulate the main existence results in Theorems l2.5l and l2.6l for the rate-dependent case 
with C — and for both rate-dependent and rate-independent cases with the positive definite mapping 
£, respectively. For the rate-dependent model with £ = 0, we assume that / satisfies either the 
coercivity condition (|1.14p or (jl.lSp . The proofs of these existence results are given in the subsequent 
sections. We note here also that the measure-valued solutions generalize naturally the notion of the 
strong solution of the problem investigated previously in [5] [51 [TO] ■ 

In Section [3| we show that for some given functions r and P the system of equations (|l.lap - (jl.ldp . 
dl.lg'P is an elliptic system of partial differential equations. Since the proof of the main existence results 
to the problem (jl.ip relies heavily on the existence theory for the equations of linear piezoelectricity, 
we use L^-existence theory for elliptic systems of partial differential equations and present the main 
properties of the solutions of the linear piezoelectricity model in full details in Section [3| 

In Section m we reduce the system p.ip to the evolution problem (|4.6p . (14. 7p . In Section [5] we use the 
Rothe time-discretization method to construct an approximating problem (|5.ip . (|5.2p and show that 
it has a unique solution. In the following Sections [Gj [7| we show the convergence of the approximating 
sequence and prove the main existence results. 

2 Statement of main results 

In this section we introduce the notion of the measure-valued solutions of the problem (jl.ll) and then 
state the main results of the work. For completeness, we give the definition of the strong solutions of 
the problem p.ip . 

Definition 2.1 (Strong solution). A function (u,(/), r, P) such that 

{u, (j)) e VF^'P' (0, T; Wq-'p' (n, X R)), (r, P) G W^^p' (0, T; L^* {n, x R^)) 

is called the strong solution of the initial boundary value problem (jl.ip . if for every t G [0,T] the 
function (u(i), 4>{t)) is the weak solution of the boundary value problem (|l.lap - (|l.ldp . ( |l.lgp with the 
given r{t) and P{t) and the evolution problem (|l.lep . (jl.lfp is satisfied pointwise. 

Next, we define the notion of the measure- valued solutions for the initial boundary value problem 
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Definition 2.2 (Measure-valued solution). A function {u,ip,r, P,t) such that 

u G W^'P' (0, T; Wo^' {n, R^)), € T4^1'P' (0, T; Wo'"' {n, R)), 
z = {r,P) e M^i'P'(0,r;i?''(f],53 x R^)), r e L^(17t,X(53 x R^)) 

with 

is called the measure-valued solution of the initial boundary value problem (ll.ip . if for every t G [0,T] 
the function (^(t), is the weak solution of the boundary value problem (|l.lap - (jl.ld|) . p.lgp with 
the given r{t) and P{t), the initial conditions (ll.lfj) are satisfied pointwise and the following inequality 

with J^it,x) = /^3xE3 (v;/)(0 nAdO holds for a.e. t e (0,r). 

Remark 2.3. We note that the integrability of the function $(s,x) = (^(p^), (_E-Vp/^(r p*)) ~ '^(p)) 
not required in Definition 12.21 We require the existence of the double integral J^^ $(s, x)dsdx, only. 

Remark 2.4. As it is discussed in AppendixIXl (|A.3|) . if we could show J^{t, x) — (r,p)f{r{t, x),P{t, x)) 
in (|2.1[) . i.e. r = <5(r,p), then the measure- valued solution (u, </>, r, P, r) became the strong solution. The 
integrability of $(s, a;) follows then automatically from the inequality (jA.ip . 

Next, we state the main results of this work. 

Theorem 2.5. Let the functions b e W^'P{0,T; LP{n,R^)), q e {0, T; LP {n,R)) and (r("),p(o)) e 
L^(r2, X M^) 6e given. Suppose that C — 0. Assume that the function g : x —i- R is convex, 
l.s.c. and satisfies growth conditions U.lffl) . Let the function f : x R^ ^ R be convex and 

such that f G C^{domf) and satisfies either the coercivity condition ^^J^ if f depends on P G R'^ 
only or U.15\) if f depends on both r £ and P £ R"^ . 

Then there exists a measure-valued solution {u, (f>, r, P, r) of the problem il.l]} . Additionally, if f 

satisfies U.15\} . then (^(p^), (^)^ integrable over fit for a.e. t £ (0,T). 
If £ is positive definite the following result holds. 

Theorem 2.6. Let the functions b £ W^'P{0,T; LP{n,R^)), q £ VV^'P {0, T; LP {n,R)) _and (r(o),P(o)) £ 
L^{fl,S^ X M"^) be given. Assume that C is positive definite and g : S'^ x R^ ^ R is convex, l.s.c. 
function, which satisfies either the growth conditions U.10\) . il.ll]) or has the form lll.9\} . Let the 
function f : x R'^ ^ R be convex and satisfy f £ (dom /) . 

Then there exists a measure-valued solution (u, (p, r, P, r) of the problem lil.l]) . Moreover, the function 

((pt)' (^) ^''"(p)) integrable over fit for a.e. t£ (0, T). 

At the end of this section we present the conditions which guarantee that a measure- valued solution of 
the problem is the strong one. The next remark motivates the introduction of the measure- valued 
solutions. 

Remark 2.7. We note that in order to guarantee that the measure- valued solution is strong one has 
to show that the inequality (|2.1[) is satisfied with F{t,x) — V (^r.P)f{f{t,x),P{t,x)). Indeed, if the 
inequality 

/iM(^:v;r.'.)-<;))--/i''(C)--a))-- 
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holds for a function g satisfying (|1.10p . (II. lip for a.e. t e (0,r), then the equivalence (jA.3[) yields 
that the measure- valued solution (u, r, P, r) is strong. For the case C — i) considered in Theorem [53] 
we just set £ = in (|2.2p . For the function g given by (|1.9p we need the following condition. In 
Theorem 12.61 we have proved that (^) — £(p) — F £ K, what yields that 5 ((^) — £(p) — -7^) = for 
the function g defined by (ll.9p . Then the verification of the following condition 

implies that the second term on the left side of the inequality (|2.ip is equal to zero, and therefore it 
ensures that the measure- valued solution (u, </>, r, P, r) is strong. 



3 Existence for linear piezoelectric models 

In section S] we reduce the system of equations (|l.ip to a single evolution equation for the vector- 
function (r, P) . This equation is just the equation (jl.lep with the functions a and E expressed through 
the functions r and P. In this section we establish the relation between the functions {a, E) and (r, P) 
using the equations (jl.lap - ljl.ldp with the homogeneous Dirichlet boundary conditions for the functions 
u and (f>. Let us suppose first that C, e and e are measurable bounded functions of x e f2. For simplicity 
we drop the time dependence of the given and the unknown function in this section. 

We use notations from section [T] and rewrite the system of equations (|l.lap - (|l.ldp as follows 

-e eiu) + eV(p J \ -q J \-er-P J 

Next, we introduce the following notations: U = {u, 0)-^, z = {r, P)^ , B = {—b, — q)^. 

Now we use the symmetric properties of the tensor C and rewrite the system of equations (|l.lap - (jl.ldp . 
p.lgP as follows 

DhA'l![(x)DkU^{x) = B,{x) + Dh£^^{x)zl{x), x e n, (3.3) 
U{x) =0, xe dn. (3.4) 

Since the entries of the mappings C, e and e are bounded measurable functions we can suppose that 
the same holds for the entries A'^^{x), x £ fl of the mapping ^ : S*^ x — S'^ x R'^. And since C 
and e are symmetric and positive definite uniformly with respect to x € fl and the terms containing e 
cancel each other in the expression A^^^ {x)-q^ -q^ with 77 e 5^ x M'^ we obtain that there exists a constant 
Co > such that the following ellipticity condition 

A^f{x)r^rll > cov^rj^ (3.5) 

holds for every 77 € 5'"^ x R"^ uniformly with respect to x G ft. 

Next, we show that the system p.3p - p.4p has a unique weak solution U € IVg^'^(57,R^ x R) for 
every given z e LP{n, x R^) and B g W^^'P{fl, x R) for 1 < p < 00. For this purpose we use the 
existence results for elliptic systems of partial differential equations. We make different assumptions 
on the entries oi A for p = 2 and p 2. If p = 2, we can apply the Lax-Milgram result for the bilinear 
form A{U,V) = J^{A^^{x)DjU^{x),DhV\x))dx to the problem ([33]) - (HH). In this case it is enough 
to suppose that p.Sp is satisfied a.e. a; G f2 and that the entries of A are measurable bounded functions. 

We prove now that for every given z e i^(f^, x R'^), B € VF~^'^(17,R^ x R) there exists a unique 
weak solution U G Wq' {n,R^ x R), which means that U satisfies 

A{U,V)^l{V) (3.6a) 

for all V = (w,V') G W^'^{n,M.^ x M), where 

A{U,V) = {A DU,DV) = {A{e{u),\'<pf,{e{v),\'^pf), (3.6b) 
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and 



1{V) = {B, V) + {£z, DV) = {b, v) + {q, V> - {£z, {siv),Vi;f). (3.6c) 

The function A : Wo''^{n,R^ x R) x Wo''^{n,M.^ x M) -> M in (IX6b| is a bilinear form. Taking 
U — V and using the elhpticity condition p.5p together with the inequahties of Korn and Poincare 
we obtain that A{U,U) > ci||C/||f 2- Since the entries of A are bounded functions we obtain that 
^)| < C2||f7||i,2||V^||i,2 and \l'(y)\ < C3||V^||i,2. Therefore, the assumptions of the Lax-Milgram 
theorem are satisfied and there exists a unique weak solution U G WQ'^{fl, x M) of the problem p.3p 
- (I3.4p . Thus, we have proved the following existence result for p — 2: 

Theorem 3.1. Let D, C be an open bounded set with dfl € C^, {'''t P) € L^{n,S^ x R'^) and 
{b,q) € VF-i'2(rj,R3 X R). Let the entries of the mappings C : S^, e : ^ e : 

be bounded measurable functions. We suppose that for a.e. x Cz the mappings C and e are linear, 
symmetric, they are positive definite uniformly with respect to a.e. x G ft, whereas e : S*"^ — > R"^ is just 
a linear mapping for a.e. x € il. 

Then there exists a unique function {u,(f>) £ WQ''^{il,MP x R), which satisfies the equations il.la]) - 
\l.ld\) with the homogeneous Dirichlet boundary conditions for arbitrary fixed t G [0,r) such that the 
estimate 

\\u\\i^2 + ||0||i,2 < c(||r||2 + ||P||2 + ||&||-i,2 + lkll-1.2) (3.7) 

holds for some constant c> 0, which is independent of r,P,b and q. 

For the case p 7^ 2 we can only prove the existence of the weak solution under the assumption that 
the functions A'^j^ix) are continuous for all a; G fi. Let 2 < p < 00. We denote fi{x) = Bi{x) + 
Dh£i-'{x)z{{x). Suppose that / G W-^^p{VL,M? x R). For every / G W-'^''P{n,B? x R) one can find 
F G LP{il, S*^ xR^) such that F satisfies fi = D^Fl in the sense of distributions and the estimate ||i^||p < 
c||/||_i.p holds. We suppose that C,e and e satisfy the assumptions of Theorem 13.21 Then since the 
assumptions of Theorem 13 . 1 1 are also satisfied, there exists a unique weak solution U G VFQ^'^(f2, R'^ x R) 
of the problem p.3p . (|3.4I1 . It is easy to prove that if C,e and e satisfy assumptions of Theorem 13.21 
then the functions A'^j^{x) are continuous for all x £ fl and satisfy the Legendre-Hadamard condition 

A'l,H^)v^VjUk>comf (3.8) 

for some Cq > and for every a: G 51, ^ G R"' and 7? G R'^ uniformly with respect to x. It is shown in O 
p. 373] that the function U belongs then to Il/Q^'^(ri, R"^ x R) and the estimate 

\\U\\i,p^n < c\\F\\p^n (3.9) 

holds with some c> independent of F. If we suppose that z G LP{n, S'^ x R^), S G W-^^P{n, R^ x R), 
then we obtain that F G LP{n,S^ x R^) and ||F|lp,n < c{\\z\\p^n + ||S||_i,p,a). 

To prove that the conclusion of Theorem 13.11 holds for 1 < p < 2 as well we use the following duality 
arguments. Let 2 < p < 00 and p* : ^ + ^ = 1. In the same way as above we prove that for any 

function / G W-^-P{n, R^ x R) there is a unique solution U G W^-^in) of the problem D^A'^^DkU^ = f, 
with the operator A replaced by A'^ such that the inequality ||?7||i,p.n < c||/|l_i^j,^o holds. Therefore 
we can define the linear bounded operator T : W~^'P{il) — )■ WQ'^{ft) by Tf — U. Then there exists a 
unique operator T* : W-'^-p' {n) ^ W^ ^' (il), such that 

[Tx,y] = [x,T*y] (3.10) 

holds for ah x G W-^-P{n) and y G W~^-P'{n). This proves that for every g G W'^-P^n) the function 
V G Wo''''{n) defined hy V = T*g satisfies DhA'^^f'DkV^ = g, in the weak sense. The uniqueness 
follows immediately. Indeed, we take arbitrary g G W~^'P (il) and U G WQ^'^(r2). The function 
/,; = DhAjiDkU^ belongs to W^^'P{il) and since U is the unique weak solution of DhAflDkU^ = fi 
we have also Tf = U. We obtain 

[U,g] = [Tf,g] = [f,T*g] ^[f,V] = - f {DhU\x),A':^{x)DuV^{x)) dx (3.11) 

Jn 
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for every U € WQ-^{fl). Since U is chosen arbitrary we get that V G satisfies DhA'^fOkV^ = gi 

in the weak sense and it is unique. From the relation T*g — V we obtain the estimate ||V^||i.p»,f2 < 
c||<?||-i.p*.r2- And since for every h ^ the estimate |jdiv /i||_-i_p. < \\h\\p* is satisfied we obtain that 
the following existence result holds for all 1 < p < oo: 

Theorem 3.2. Let Q C M!^ be an open bounded set with dn G C\ (r, P) G LP{n,S^ x R^) and 
ib,q) e VF-i'P(ri,]R3 X M) wi^/i 1 < p < oo. Lei t/ie entries of the mappings C : , e : ^ 

IR'^, e : ^ be continuous functions of x ^Q,. We suppose that for every a; G O the mappings C 
and e are linear, symmetric, they are positive definite uniformly with respect to x € ^, and the mapping 
e : ^ M.'^ is linear. 

Then there exists a unique weak solution {u,(j)) G W^'P{fl,M.^ x M) of the problem /il.la]} - [l.ld\} with 
the homogeneous Dirichlet boundary conditions for arbitrary fixed t € [0, oo) and the estimate 

lkl|i,P + UKp < c(||r||p + ||P||p + + ||g||_i,p) (3.12) 

holds with a constant c> 0, which is independent ofr,P,b and q. 



4 Reduction to the evolution equation 

In this section we show that the function (cr, E) can be expressed conveniently through the function 
z — (r, P) in such a way that after substituting (cr, E) into the equations (|l.lep and (jl.lfl) the problem 
p.ip is reduced to an evolution problem for the function z. 

Let us suppose that the function {r,P) is known and belongs to LP{VI,S^ x R'^), 1 < p < oo. We 
consider the equations (|l.lap - (ll.ld|) . p.lg[ ) and suppose that the entries of the mappings C : 
S^, e : M.^ —i' M.^, e : M."^ are continuous functions of x G fl. Since the assumptions of Theorem 

13.21 are satisfied we get that for every given z = {r,P) G LP{il,S'^ x M.^) this problem has a unique 
solution U = {u,(j)) G Wl"^{Vt,S^ x R'^). Let us decompose U = Ui^ + Ub, where Uq — (uo,0o) is 
a solution of the problem (|l.la|) - (jl.ld[) with {b,q) = and {r,P) ^ 0, and IJb — {ub,4'b) satisfies 
(frTa|) - (frTd| with (&, q)^0 and (r, P) ^ 0. 

It follows from Theorem 13.21 that for all 1 < p < oo the following estimate holds for the functions uq 
and (j)o 

||"o||i,p + ||0o||i,p<c(||r||p + ||F||p). (4.1) 

Next, we define a linear operator Q : LP{n, x R^) -> LP{n, x R^) by 

Q{r,Pf = {e{uo),Dof, (4.2) 

which is bounded due to (|4.1|) . It turns out that Q is a projection operator. To this end, we consider 
the functions Eq = Eq — r and Dq — Dq — P and rewrite the equations (|l.lcp . (jl.ldp as follows 

{ao^Eof ^ViEo,Dof, (4.3) 
where I? : S''^ x R'^ -> S*^ x M'^ is the operator defined by 

It was shown in [5| that V is symmetric and positive definite. 

L/emma 4.1. Let the vector {r, P)^ G LP{n, x R^) be given. We define the linear mapping Q = Qp : 
LP{n,S^ X R3) ^ LP{n,S^ X R3) by where {e{uo),Dq)^ .satisfies the problem fTTT^j - flT^) with 

{b,q) = 0. Then Qp is the projection operator, which is adjoint to the operator Qp- with respect to the 
bilinear form [zi,Z2]v = {T^zi, Z2)p,p* ■ 

Proof: The operator Qp maps the elements of the space LP{Q,,S^ x R^) into the subspace Ti. = 
{w = {E{uf)),DoY ■■ uo G Wo'''(r2,R3), Do G LP{n,R^) : div Dq = 0} of the space LP{n,S^ x R^). 
Since (j4.3p contains only the differences e{uq) — r, Dq — P we obtain from the uniqueness of the 
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solution of the problem (|l.lap - (jl.ld|) with homogeneous Dirichlet boundary conditions that for every 
w e LP{fl, X M'^) the operator Qp is the projection operator, namely, Q^w — QpW. 

Now we show that Qp- is the adjoint operator to the operator Qp with respect to the bilinear form 
[•, ■]v, i.e. we show that for every z = (r, P)^ g LP{n,S^ x R^) and z* = {r\P*Y 6 L^' [VL.S^ x E^) 
the equality 

[Qp*z*,z\t) = [z*,Qpz\-D (4.5) 

holds. 

Let us denote the images of Qp*z* and QpZ as Qp*z* = (£q,Z)q)-^ and QpZ = {eQ,Do)'^. We have 
{ao,Eof = V{{sa,Daf - (r, P)^) = ViQpZ - z) 

and in the same way {aQ^E^y = 'D{Qz* — z*). We show now that [{Qp* — I)z* ,Qpz]v — and 
[{Qp ~ Qp*z*]v — are satisfied, it would imply (|4.5p . We obtain 

[{Qr^I)z*,Qpz]v = {V{Qp,z* ^z*),Qpz) = {{a;,E*of,{eo,Dof) = {a*,e{uo)) 
+ {-V(f,lDo) = - (div a*o ,uo) + {(f>*Q, divDo) = 0. 

In the same way one can show that [{Qp — I)z, Qp-'Z*]xi — holds. Then we have proved that Qp is the 
projection operator, which is adjoint to the operator Qp- with respect to the bilinear form [■,■]©. □ 

Now, let us define M V{I - Qp), z := {aB^Es)'^, z := {r,PY and z° := (r°,P°)^. Inserting 
the expression (ctoj-Eo)"^ = 'D{Qp — I){r,P)'^ into the equation (jl.lep with the initial conditions (jl.lf|) 
yields that equations (jl.lej) - ()l.lf|) can be rewritten in the following abstract form 

ZtedIg{-Mz^Cz~dIf{z)+z), (4.6) 
^(0) = 2°, (4.7) 

where /g, // : L^(r2, 5^ x R^) ^- R are functional defined by (|A.7p . 



5 Existence and uniqueness for a time-discretized problem 

We show the existence of measure-valued solutions using the Rothe method (a time-discretization 
method). In order to introduce a time-discretized problem, let us fix any to € N and set h := From 
the assumptions for the functions b and q we can conclude that z e LP{n, x R^). We set 



nh 

3 w T03\ 



z{s)ds e LP{n,S-' xW), n==l,...,2'^ 

(ri-l)/i 



Then we are looking for functions z"^ e L^{^, x R'^) solving the following problem 

~n ^n— 1 

" edi,{j:-J, (5.1) 

= (5.2) 

with 

S;^, -M^zZ - dlf{zl,) + e ^3 ^ j^3)_ (5 3) 

where 

:= (2?(J ^Q2)+C+ —I) : L^{n, x R^) ^ L^^-^^ 5.3 ^ ]^3)^ 

TO 

To show that the discretized problem has a solution we need Mm to be positive definite. This holds due 
to the term ^/ even if C is only positive semidefinite. Therefore we consider here rate-dependent case 
with C = and rate-independent case with £ > simultaneously and suppose that C is only positive 
semidefinite. Recall that the functionals // and Ig are proper, convex and lower semi-continuous (see 
Section]^. We want to show that the equation (|5.1I) can be rewritten as 

d^{zl) 9 zl, (5.4) 
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where 

The functional ^\\M^J'\\\l : L^{n, x M^) ^ ^ is proper, convex and lower semi-continuous. 

Indeed, since Mm is bounded and positive definite operator, then it is maximal monotone by Theorem 
II. 1.3 in [T]. Since the operator Mm is also self-adjoint, one has that Mm = by Proposition II. 2. 7 
in [1]. All other properties of $ follow from its definition. The last thing which we have to verify is 
whether the following relation 

9* = dig-. + 9$ + dif 

holds. By the definition of we conclude that the domain of $ is equal to the whole space x 
M^). By condition (II. lip the domain of the functional /g* is also the whole space 
the rate-dependent case. In the rate-independent case the domain of /g* also coincides with L^{ft, x 
M.^). Therefore, condition (jA.6| is fulfilled and, since all functionals are proper, convex and lower 
semi-continuous, Proposition IA.3I gives the desired result. With the relation (|A.2p in hands the last 
observation implies that 

dom(a^') = dom(9//). 

Since $ is coercive in L'^{^, x K^), which obviously yields the coercivity of the operator A = d'i' 
is surjective by Theorem IA.21 Thus, we conclude that for every fixed m S N and n — 1,...,2™ the 
problem I^J^f, (H^]) has a solution zl'^ e L'^{n,S^ x R^) for every given G LP{n,S^ x M^) and 
G L'^{^l,S^ X E'^). The solution zj^ is also unique. Indeed, suppose there are two functions zi and 
Z2, which satisfy the equation (|5.4p for a given Zm- We substitute the functions zi and Z2 into (|5.4p 
and consider the difference of both equations. Then using the monotonicity of dig* and dIf we obtain 
that 



(M„(zi - Z2), zi - Z2} < 0, 



which together with the positive definity of Mm implies that the solutions coincide. 

Rothe approximation functions: For any family {fm}n=o,...,2™ of functions in a reflexive Banach 
space X, we define the piecewise ajfine interpolant G C{[0,T], X) by 

^mit) := in- 1)) C + - {n-l)h<t< nh (5.5) 

and the piecewise constant interpolant G L°°(0,T; X) by 

W := C for {n - l)h <t<nh, n=l, 2", and ^^(O) := fm- (5-6) 
For the further analysis we recall the following property of ^m and ^m'- 

||Cm||L=(0,T;X) < ||?m||L=(-/!,T;X) < (^/l||Cmllx + II Cm 1 1 (0,T; Jf) j ) (5-7) 

where is formally extended to t < by and 1 < s < cx) (see [13]). 



6 A-priori estimates. 

Rate-dependent case. We suppose that g satisfies the conditions (|1.10p and (jl.lip . Let us fix 
m G N and n = 1,...,2'". Since the problem (|5.ip . (|5.2p has a unique solution, we obtain with the 
Young- Fenchel property (see Appendix 

yTL ..n—\ , n _n— 1 , 

J ( m m \ I T {^n \ / m m y^n \ 



11 



which together with the relation (|5.3p impUes 

•^m ~ ^rrT^ 1 / 1 \ 

+ ^(C - z':;^\dif{zl)) = (f^i^, z;;). (6.1) 

We note that {z!^ — z^"^, > (t>{z^) — (j){z^^) holds for any convex functional 4>. Therefore, 

multiplying (|6.ip by /i and summing the obtained relation for n — l,...,l for any fixed / S {1, ...,2™} 
we derive the following inequality 

n— 1 n— 1 

< ^ (ii (M + /:)i/2,o||2 + i.||,o||2) ^ ^^(^0) ^ 1^ ^ II n,, ||i;j,||,. (6.2) 

n=l 

Applying the conditions (jl.lOp and (|l.lip to the terms, which contain Ig and Ig, and the Young 
inequality with e < di to the last term in (|6.2p we obtain 



Hd, -e)J2 f^'"' H^: + /^ci ^ ||s;^||^ + - (||(M + cy/'z'jil + -1141I2) + If i^L) 

71—1 n—1 
< i(||(M + Cy^'z^Wl + -Wz^Wi) + //(z") + C^hj^UrnW; + id2 + C2)\nT\. (6.3) 



Now, taking Remark 8.15 in |13j into account and using the definition of Rothe's approximation 
functions we rewrite (16.31) as follows 



(rfi - emz^\\i:,,^+c,\\f:,ji,,^ + ^{\\iM+cy/^z„M\i + ^ji^mimi) 

+ If{zm{t)) < l{\\M'/'zY2 + ^\\zli) + If{z') + {d, + c^mrl + C,fzm\\l,n,. (6.4) 

Since — > z in LP(f7T), the last term in (|6.4p is bounded by a constant. The estimate (|6.4p implies 
that 

{z„}„ is uniformly bounded in W^^p' {0,T; L^' {n, x R^)), (6.5) 

{E„,}^^ is uniformly bounded in LP(rjT, x R^), (6.6) 

{(M + /:)i/2z„J,„ is uniformly bounded in L°°(0, T; ^^(f]^ x R^)), (6.7) 

^ z,n \ is uniformly bounded in L°°(0, T; ^^(f^^ ^3 ^ j^s^^^ j-g 



{/(2m)}m and {/(z„)},„ arc uniformly bounded in i°°(0, T; L^fi, R)). (6.9) 

We can improve the estimates (j6.6l) - (|6.8p and show that is uniformly bounded in the space 

L°°{0,T; LP{fl, S"^ X R'^)) and the sequences {{M + Cy^'^Zm}m and {-^Zm}m are uniformly bounded 

in W^''^{0,T; L^{fl, S"^ x R'^)). Indeed, multiplying (15.11) by the term ^""^f "* and integrating over fl 
we obtain 



, _n _ n—1 yn _ yn — 1 \ 1 / v 1 



= - a/g(E;:j,E:^ - Kn')> TilgiK^) - IgiKn'))- (6.10) 
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Then we use that dif is a monotone operator and estimate the left side of (|6.10p from above 

-'m m ni m \ / m m (^^f -\- d -\- — — ) 



ft, ' m h 

n _ n-1 flT^f .,n \ _ -,"-1 \ , _ _n-l _ -ri-1 

'-'-^IX'^m) "-^J\'^m )\ , / 



<^ _ ||||^/ -|- 11"^ — II 77?. m 1 1 p _j_ ^ "^m '^m '^m '^m ^ (6 11) 



Now we combine (j6.10p and (16. lip , multiply the obtained relation by h and sum it up for n — 1, 
and any fixed / € {1, 2^"-}. We obtain 

^ Tl TL — 1 1 TL Ti — "1 C 

ft^(||(Af + /:)V^ "-7- ii^ + n^ "--"" 11^)+/ g(sL)dx 

< /,(S(0)) + ft^ ( " ^"^ , " ), (6.12) 

n=l 

which imphes the estimate 

||(Af + + ||^z„.,|||^^ + < /g(E(0)) + ||z™t|ir,n,pmt||p.n,- (6-13) 

Since by (16.51) the rigth side of (|6.13p is bounded we obtain 

is uniformly bounded in L"^ {0,T; LP{n, x R^)), (6.14) 
{(Af + Cy/'^Zrnjm is Uniformly bounded in W^'^{0, T; ^^(f], 5^ x R^)), (6.15) 

1 is uniformly bounded in W^-^{0, T; L^{n, x R^)), (5.16) 



m 

In conclusion we note that if the function / depends only on P and satisfies the coercivity condition 
pm)) . we would obtain from (HH) 

{P,n}m and {P,n}m arc uniformly bounded in L°°{0, T; L^{n, x M^)) (6.17) 

and if / depends on both r and P and satisfies the coercivity condition (|1.15p , we would get then from 

and {z™}™ are uniformly bounded in L°°(0,T; LP{n, x R^)),. (6.18) 

Here we emphasize that the estimates of this paragraph hold for positive semi-definite operator C. 
Suppose now that C is positive definite. Then the estimate (j6.15l) immediately implies that 

{z™}™ is uniformly bounded in W^^^{0,T; L^{n, x R^)) (6.19) 

without coercivity assumptions for the function /. 

Rate-independent case. Now we suppose that g is defined by (|1.9I) and C is positive definite. The 
proof runs the same lines of the second part of the previous paragraph except some slight changes. We 

multiply (|5.ip again by the term — , integrate over D, and use that /g(S]^) =0, m,n e N to 

obtain 

~n _n— 1 yn yn — 1 \ 1 / v 1 

" ^"^ , " ^ " ) = /^(5/s(S;j.), S:;. - 7 > J^ilrAKn) lAKn')) = 0- (6-20) 

It follows from ((00)1 that 

1 \ ~ri „n — 1 n 1 , l , 

A/ + £ + - j - ^"^ , " ) + ^(a//(z;^j - difizi^-'), - z--^ 

I -,n— 1 in— 1 . 

^ / \ 2]^^ 
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Since M is positive semidefinite, dif is monotone and C is positive definite we obtain the estimate 



m in 



< c 



2, SI 



2, a 



which after multiplying by h and summing for n = 1 , . . . , Z implies 



(6.22) 



Since the equation (|6.1I) holds also rate-independent we use it with the same arguments as in the 
previous paragraph to obtain the following estimate 



t pt 



< i(||(M + C)'/'zY2 + + //(^") + \\z„u\knA~Zm\\ 



2,07 



(6.23) 



Since the sequence {zm}m converges strongly in W^'P{0, T; LP{il)) and the domain of Ig is the bounded 
convex set K we obtain in the rate independent case that the estimates (|6.22p and (|6.23p imply 



{z„}™ is uniformly bounded in W^'^{0,T; L'^{i}, x K^)), 

{f^rn}„^ is uniformly bounded in L°°{0,T; L°°{n, x R^)), 
{/(2m)}m is uniformly bounded in L°°{0,T; L^{n,R)). 



(6.24) 

(6.25) 
(6.26) 



7 Existence of measure-valued solutions 



Based on the results of the previous sections we are able now to prove the main existence results of this 
work, Theorem [23] and Theorem 12^ 



The proof of Theorem [2751 

Proof. In a similar way as in the beginning of section [5] we use the equality (j6.ip to derive the following 
inequality 



/ g*{zjnt{x,s))dsdx + / / g{Em{x,s))dsdx + {MZm,Zmt) 

+ 1 f{z„i{x,t))dx - 1 f{z°{x))dx< ■^\\z'^\\l + {Zjnt,Z„i)2^^^ 



2,Ot 



(7.1) 



Using a-priori estimates from section [B] we want to pass to the limit in the inequality (j7.ip . 

First, we obtain by (|6.5p that, at the expense of extracting a subsequence, the sequence {zm}?? 
converges weakly in the space M^^'P {Q,T;LP (il)) to some z. Next we claim that the sequence {zm}n 
converges weakly in (fir) and the weak limits of {zm}m and {zm}m coincide. Indeed, using (|6.5 
this can be shown as follows 



IP 



nh 



n=l"'("-l)'i 



(z" - 



i^t — nh 



dt 



p* + 1 ^ 



hp 



dZr, 



dt 



(7.2) 



which implies that Zm — Zm converges strongly to in U (fir)- Then the sequence \zm^ra converges 
weakly in L^' {^t) to the same weak limit z € W^-p {0,T;LP (Q)) as {zm}m- 
In the same way, using (|6.15p and (|6.16p . we obtain that 



M(i 



Zm) and — (z„i — z„i) — > in i (ilr)- 
m 



(7.3) 
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Since the elements of {E,„}meN belong to the space LP(ilT), we can use the definition (jA.7[) to define 
the convex functionals Ig : U'{D,t) — ?> M and Ig* : (fit) R by ig*{z) = g* {z{s, x))dsdx and 
ig{z) = 5(2(5, a;))(is(ia;, respectively. Then the functionals satisfy /g* — {Ig)* . Since z^t converges 
weakly to zt in {^t) we obtain 



liminf Ig, (zmt) > Ig' (zt). (7.4) 

m— voo 



Due to (|6.6p there exists a subsequence of which converges weakly to S in LP{ilT)- Therefore 

we get 

liminf 4(Sm) > /g(^)- 

m— ^00 

Our next goal is to compute the function S. To this end, we consider the sequence E„i — —Mzm + 
^z„i - df{zm) + Zm,_'m eN. 

By constructions, z„i converges strongly to z in U'{Q.t)- It follows from (|6.16[) that the sequence 
{-^z„i\m is uniformly bounded in L'^i^lx). By (|7.3p the sequence {-j=z„i}m is also uniformly bounded 
in IP'{Q,t) and therefore converges strongly to in L'^{^t) ■ 

Moreover, it follows from (|6.15p and (|7.3p that the sequence {Mzm}m converges weakly in L'^i flx) and 
the weak limit coincides with the weak limit of the sequence {Mzm}m — {(^'")}m- By Theorem 13. 21 we 
get that for every z^ & W^'^{0, T; L'^{fl)) there exists a solution (cr,„, E„i) € W^''^{0, T; L'^{fl; x ] ' 
of the problem 



-div ffrn = 0, div Dm = 0, (7.5) 



-V{I - Q2)z„, = -Mz„ (7.6) 



with homogeneous boundary conditions for the functions (urm(t'm)- 

Now, let us define the operator Mp* ~ T>{I — Qp*) : (il) [Vi). Due to the uniqueness of the 

solution of the problem (|7.5p - (|7.6p . instead of the operator M we can consider the operator Mp. in 
(|7.6p . We note that the operator Mp* is the extension of the operator M on L^(r2). Thus, because of the 
linearity of the operator Mp. , the sequence {(^")}m converges weakly in the space Vl^^'^(0, T; L^(r2)) 
to Mp.z = (^) , which is the solution of the problem ()7.5p . (|7.6p corresponding to the function z. Since 
the sequence {—Mzm + —Zm + Zm}m converges weakly in L'^{VLt) and {S,„}„i converges weakly in 
L^{^t): the sequence {dl f{z„i)}m converges weakly in L'^{n,T) to some T e £^(0t). By Theorem 
IB. 31 there exists a Young measure r S Lu,(r2T; A^('S'^ x R'^)) associated with the sequence {zm}m such 
that F{t,x) — Jg3^^3 df{£,)Tt^x{d(,) for a.e. {t,x) € 51^- Thus, we get that Em converges weakly to 
—Mp*z — + z in U'{VIt) and that the inequality holds 

liminf /„(E™) > iJ~Mp,z - F + z). (7.7) 

Since J{z) = f{z{x))dx is a convex functional on (ft) and Zm{t) converges weakly to z{t) in 
LP (ri), we get that liminfm->.oo Jq f{zm{x^t))dx > f{z{x,t))dx. Next, for a.e. x £ VL we have that 
z{x, ■) e iy^'P*(0,T; 5*3 X R3)_ Let us fix x £ From JlH) we conclude that the set /(z(a;, [0,r])) C 
dom/. Since / S C^(dom/) and z(a;, [0,T]) is a compact subset of dom/, we obtain that /(z(-)) : 
[0,T] — >■ R is an absolute continuous function and therefore f{z{t)) is almost everywhere strongly 
differentiable with /(z(i)) - f{z°) = J*{df{z{T),ztiT))dT for a.e. t e (0,T). It follows from (HI]) that 
J^{f{z{x, t)) — J {z^ {x))dx < 00 and hence ^^{df{z{x, r), zt(a;, T))dTdx < 00. Therefore we ge10 

limM f if {z„^{x,t))~ f{zy{x))dx> f f {df{z{x,T),zt{x,T))dTdx. (7.8) 



m— f 00 



n JnJo 



^The existence of the integral f^{df{z{x, t), zt(x, T))dxdT is an open problem. 



15 



Now, we estimate lim inf „i_5.oo {Mzm, z„it)2 Ut fr'^™ below. To this end, we note first that (|7.3p implies 
that lim inf „i_5.oo {Mzm, Zmt)2 n ^ Ihii inf „i_>.oo {MzrmZmt)2 o • The equations (j7.5l) - (j7.6p yield that 

a:) 

We introduce the functions Wm — (^"^ ) • Then in terms of Wm the last term can be rewritten as follows 

{V-'w,n,W,nt)2,n, = 1 1 ' ) ^ W 1 1 2 " 1 1 (2?"' ) ^ "'(0) 1 1 ^ 

where w{0) = Wm{0) does not depend on m since ^^(O) — z^. The relation ~ {Mzjn}jn implies 

that w„i ^ u> in W^'^{0,T; L^{il)) and since is a positive definite symmetric operator, we obtain 
that (T)^^)iwm{t) {'D~^)^w{t) in L^{il). Since w i— )■ is a proper, convex and l.s.c. 

functional on L^(fi), using Lemma \AA\ we obtain 

liminf (2?-im„.,u;™t)2j. = liminf (|| (2?-i)^u;„(i)||2 - || (2?-i)5u;(0)||2) 

> - \\{V-')iwml = {'D-'w,w,)2,^^. (7.10) 

Let {u,(j),CF,D) be a solution of the problem (|7.5p . (|7.6I) with homogeneous boundary conditions corre- 
sponding to z. And let us assume first that / satisfies the coercivity condition (jl.lSp . Then we have 
that Zm^ z'm LP{VtT). Combining ([7?^)) . (|71II| we obtain 



m—^QO 



JnKyE{V<t>)J'\ A Jjd^d.s + {Mz,z,)^^^,^^^ (7.11) 

Since z € L°°(0, T; LP(17)), then also (£;(v^)) e L°°(0, T; LP(fi)) and we obtain that the first term on 
the right side of CHI {{Eivc^)),{'%:'^))p,p^,n^ = 0. Therefore 

Iminf {Mz^,z^t)2^nt ^ (^'^^' ^*)p,p*,n* ■ (7-12) 

Altogether we obtain 

ig,{zt)+ig{-MpZ-T + z)<{zt,-MpZ + z)pp,^^^- J J^{zt,df{z))dtdx. (7.13) 

Suppose now that f{r,P) — f{P) and satisfies the coercivity condition (|1.14p . The coercivity con- 
dition implies that P e L°°(0, T; L^(i7)). Since df{z) = (gj^p-j) and the sequence (!]„)„ is uni- 
formly bounded in Lf{nT), we obtain that a G LP{ilT) and {a.et)^^^ — 0. Now since E,D — P e 
W^'^{0,T; L^{il)) we integrate by parts the term {E,Dt — Pt)^^ with respect to t to obtain 

{E, A - Pt)^,, = {E,D- P)J, - {EuD - P)^^. 

Since P e L'^{n) and £» - P e L'^{n) for every t e [0, T], then also D e L'^{n). Therefore we obtain for 
every t e [0,T] that {E,D)^-^ = 0. And since Et = E{V<l)t) G L'^{Vt), we obtain for a.e. t £ [0,T] that 
also {Et, D)^ =0. Thus we obtain that 

(£;, A - Pt)n, = -{E,P)n\l + {EuP}n,- i^.U) 

Now we show that the function x — > Jq{E, Pt) G L^{Vl). Indeed, we have that for every fixed x G 17 the 
function [E, P) is absolutely continuous and 

d 



{E,P)\l^ j^{E,P)ds^ j^{EuP) + {E,Pt)ds 
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It follows from the estimate (|7.14p that J^{E, Pt)dsdx exists. Therefore we get that the equality 
In /o ^* ~ Pt)dsdx = — Jq{E, Pt)dsdx is satisfied and with to oo the estimate (|7.ip takes the 



ig,{zt)+ig{-Mp,z-T + z)<{zt,z) f f (zt,Mp.z + df{z))dtdx. (7.15) 

' ' Jn Jo 



form 



This completes the proof of Theorem [23] □ 
Now we prove the existence result for the case when £ is a positive definite operator. 

The proof of Theorem [2761 

Proof. We start again from the inequality 

g*{z„uix,s))dsdx + / g{f]jrt{x, s))dsdx + {{M + C)zm: z^t)2 
Jn Jo Jn 

+ I f{Zra{x,t))dx- I f[z\x))dx < ^\\zY2 + {Zmuln)2,n,- (7-16) 

Using (|6.19p for the rate-dependent and (|6.24p for the rate-independent case we obtain based on the 
same arguments as in the \1 .2\ that up to a subsequence the sequence {zm — Zm,}m converges strongly 
to in L'^{ilT) and therefore the sequences {zm}m and {z™}™ converge weakly to the same limit 
z e PF^'^(0, T; L^(r2)) in L'^{D,t) and {zmtjm converges weakly to zt e i^(r2T). 

The functionals Ig : L^{flt) — > M and Ig* : L^iflt) — ^ R defined by Ig* (z) = /J /^^ g*{z{s, x))dsdx and 
ig{z) = jQg{z(s,x))dsdx, respectively, are convex, proper and l.s.c. for both rate-dependent and 
rate-independent choices of the function g. Obviously, it holds Ig* — (Ig)* ■ Since Zmt converges weakly 
to Zt in L'^lflx) we obtain 

liminf ig. (zmt) > Ig' (zt). (7.17) 

Due to (|6.6p and (|6.25p there is a subsequence of {T,m}m, which converges weakly to E = — (M +C)z— 
J^-l-z in L'^{Qt) with the function J" G L^(17t) such that J"(t, x) = /gsxH^ ^fiO'^t,x{d£,) for a.e. {t, x) £ 
Or with the Young measure r £ L^(yiT]M.{S^ x R"^)) associated to the sequence {zjn}m- Then we 
obtain 

liminf /„(E™) > ig{-{M + C)z-J- + z). (7.18) 

Since Zm(t) converges weakly to z(t) in L^(r2), then liminfm_j.oo f{zm{x,t))dx > J^^ f{z{x,t))dx. 
Let us fix xe n. For every fixed X G we have z(x, •) £ W^'^{0,T; xM.^). Since / G C^S'^ x R^, R), 
we obtain that /(z(-)) ; [0,T] ^ R is an absolute continuous function and therefore f{z{t)) is for a.e. 
t £ (0,T) strongly differentiable with f{z{t)) - f{z°) = J^{df{z{T),zt{T))dT. It follows from (dH) and 
(|6.26p that j^{f{z{x,t)) — f{z'^{x))dx < oo and hence J^-^ J^{df{z{x,T), zt{x,T))dTdx < oo. Therefore 
we have 

limM f if {zm{x,t))~ f{z"{x))dx> f f {df{z{x,T),zt{x,T))dTdx. (7.19) 



m—^oo 



n Jo 

2/ 



Due to the weak convergence of the sequences {zm)m and {zm)m in the space L (Qt) we obtain 
liminf ((M + £)z„, z^t), o, = liminf ((A/ -I- /:)z™, z^t), o, = liminf(||(Af -|- /:)^Zm(i)||^ 



m— ^oo 

2 



- \\{M+c)^zmi) > \\{M+c)-2zmi - \\{M+c)-^zmi = m+c)z,zt)2,,, (7.20) 

Altogether we take the limit inferior on the left side and the limit on the right side of (|7.16p and obtain 
with (|7T7| - (|7?^ the following inequahty 

ig.izt) + Ig{-{M + C)z-T + z)<{zt,~{M + C)z + z)^,.- f ( {zt,df{z))dtdx. (7.21) 

Jn Jo 

This completes the proof of Theorem 12.61 □ 
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Remark 7.1. In the proof of Theorem [53] and Theorem l2.6l the function (r, P) is obtained as the hmit 
of the weakly convergent sequence {rn,Pn) and the measure r is the Young measure associated with 
this sequence (r„, P„) (see Appendix |B] for the main properties of Young measures, if needed). As we 
have mentioned in Remark 12.41 the measure- valued solution (u, (j),r, P) is strong if = V(r p)/(r, P) 
in (|2.ip . This holds if the mapping ^(^r,p)f is affine. However, this case has no practical applications 
and therefore beyond our interests. 



A Convex analysis 

In this section we briefly recall some basic facts about convex functions, their subdifferentials and the 
surjectivity results for them. 

Let be a reflexive Banach space with the norm || • || , be its dual space with the norm || • ||*. The 
brackets (•,•) denote the duality pairing between V and V* . By V we shall always mean a reflexive 
Banach space throughout this section. 

For a function 4> : V M. the sets 

dom((/)) ^{v eV \ (t>{v) < oo}, epi((/>) ^ {{v,t) eV xR \ (f){v) < t} 

are called the effective domain and the epigraph of 0, respectively. One says that the function is 
proper if dom((/)) ^ and (f>{v) > ~oa for every v V. The epigraph is a non-empty closed convex set 
iff (/) is a proper lower semi-continuous convex function or, equivalently. iff is a proper weakly lower 
semi-continuous convex function (see [16, Theorem 2.2.1], if needed). 

The Legendre-Fenchel conjugate of a proper convex lower semi-continuous function : y — > R is the 
function (/>* defined for each v* e V* by 

(j)*{v*) = sup{{v*,v) - 0(w)}. 

The Legendre-Fenchel conjugate (j)* is convex, lower semi-continuous and proper on the dual space V* . 
Moreover, the Young-Fenchel inequality holds 

Vv e V, Vv* e V* : (f>*{v*) + <l){v)>{v*,v), (A.l) 

and the inequality (j) < ip implies ■0* < 0* for any two proper convex lower semi-continuous functions 
■0, : — > K (see [16", Theorem 2.3.1]). Due to Proposition II. 2. 5 in [1] a proper convex lower semi- 
continuous function satisfies the following identity 

int dom(0) int dom(50), (A. 2) 

where 90 : F — > 2^ denotes the subdifferential of the function 0. We note that the equality in (jA.l|) 
holds iff V* e 90(w), i.e. together with (|A.ll) 

V* e a0(w) ^ 0* [v* ) -t- 0(w) < {v* , u) , Vu e Vu* € V* . (A.3) 

Remark A.l. We recall that the subdifferential of a lower semi-continuous proper and convex function 
is maximal monoton^ (see [Ij Theorem II. 2.1]). 

The next surjectivity result on subdifferentials of convex functions is one of the key tools in the proof 
of our main existence result. 

Theorem A. 2. Let A := d<j), where (j) : V ^ ^ is a proper convex lower semi- continuous function on 
V . Then the following conditions are equivalent 

hm M=oo; (A.4) 
R{A) = V* and A-^ is bounded. (A.5) 



■^A monotone mapping A : V 2^ is called maximal monotone iff the inequality 

{v* - u* ,v - u) > V n* G 

implies v* G A{v). 



18 



Proof. See PQ Theorem II. 2. 6], for example. □ 
To state our next resuh, we recah that the relation 

del) + dip ^ d{4> + V) 

holds for any two convex functions ip and (j), if there exists a point in dom((/)) n dom('0) where (j) is 
continuous (see |15l Proposition II. 7. 7]). Then, since a proper convex lower semi-continuous function is 
continuous on the interior of its domain (p] Proposition II. 2. 2]), we get the following important result. 

Proposition A. 3. Let 4> be a proper convex lower semi- continuous function and ip be convex. Suppose 
that 

int dom((^) n dom{i;) ^ 0. (A.6) 

Then 

We will use the following chain rule 

Lemma A. 4. Let (f> : H ^ Moo be proper, convex and l.s.c. on H. If z,Zt G L^(0,T; H) and if there 
exists g G L^{0,T; H) with g G d(f>{z) a.e. on [0,r], then 4>iz) is absolutely continuous on [0,T] and 



holds for every h G d(f>(u) a.e. on [0,r]. 



Convex integrands. For a proper convex lower semi-continuous function </) : R'^ — > M we define a 
functional on Lp(17,M'=) by 

r . ^ f Jn<l^ivi^))dx, m G Ll(f7,R'=) 
I +CXD, otherwise 

where is a bounded domain in with some G N. Due to Proposition II. 8.1 in [T3] the functional 
Iip is proper, convex, lower semi-continuous, and v* G dl^{v) iff 

V* e LP' veLP{n,R'') and v* (x) e d(t){vix)), a.e. 

Due to the result of Rockafellar in [TTJ Theorem 2] the Legendre-Fenchel conjugate of is equal to 
I4,' , i-e. 

where </)* is the Legendre-Fenchel conjugate of (j). 



B Young measures 

Let E C M" be a Lebesgue measurable set with ^{E) < 00. We denote Co(M'*) = 
Let 7W(R'') of signed Radon measures with bounded total variation. There is a one-to-one corre- 
spondence between the dual space of Co(K'^) and the space AA{W^), such that v G A^(R'^) defines a 
linear continuous functional on Co(M'') in the following sense {v,f) = f{x)v{dx) for / G Co(M'^) 
and ||i^||7V((R'') — '^uP||/||oo<i K*^'/)!- C)ne says that v G A^(M'*) belongs to Prob(IR'') \i v \s & prob- 
ability measure. The mapping t : E A4{M.'^) is said to belong to the space L'^{E,M{M.'^)), 
if for all / G L^{E;Co{R'^)) the function x {r^^f^x,-)) = J^^ f{x, X)dTx{X) is measurable and 
\\'^\\l^{E;M(«'^)) = ess sup^g£;||r2,||yn(Kd) < 00. 

Definition B.l. A Young measure t : E A^(M'') is an element of the space L'^{E,M{R''')), such 
that G Prob(M'*) for ^-a.e. x G E. 
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Theorem B.2. Let {m„}„ : E ^ R'^ be a norm bounded sequence in L^{E;R'^). Then there exists 
a subsequence {wntlfc of the sequence {w„}„ and a Young measure r G L^(S, A^(M'^)) such that the 
sequence {t„j.}a; konverges to r in L^{E,M.{R'^)). 

Theorem B.3. Let : -E — > R'^ be a sequence of measurable functions and {t„}„ be a sequence 

of Young measures associated to functions {un}n such that Tn ^ t. Let $ : M'' — > M 6e continuous and 
suppose that {^{un)}n is uniformly integrable. 
Then 

J E JR'i 

and converges {L^,L°°) to w, where 

w{x) = / '^{£,)Tx{d£,) for ^ - a.e. x G E. 
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C Notations 

(•,•), I • |: skalar product and norm in or 

(•, •), II • II2: skalar product and norm in Z/^(0) or Z/^(f]) 

(•, ■)p^p*, II • llpi bilinear form on L''(f7) x L^ and norm in LP{fl) 

(''')fe,p,n' II ■ IU,p,f2- skalar product and norm in W'''P{n) 

[•, •]: is a bilinear form on Wo'^'in) x W-^'P' (n). 
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